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TIME SPLITTING SCHEMES AND MEASURE SOURCE TERMS 
FOR A QUASILINEAR RELAXING SYSTEM 



Several singular limits are investigated in the context of a 2 X 2 system arising for instance 
in the modeling of chromatographic processes. In particular, we focus on the case where 
the relaxation term and a L 2 projection operator are concentrated on a discrete lattice 
by means of Dirac measures. This formulation allows to study more easily some time- 
splitting numerical schemes. 

Keywords: relaxation schemes, nonconservative products, conservation laws. 

1. Introduction 

Relaxation is a phenomenon which appears in a wide variety of physical situa- 
tions. In gas dynamics, it occurs when the gas isn't in thermodynamic equilibrium. 
In elasticity, it is usually referred to as a fading memory mechanism. Hyperbolic 
conservation laws with relaxation are also served as discrete kinetic models. Re- 
laxation approximations to scalar conservation laws can be considered for purely 
numerical purposes too; general references are e.g. 9 > n . 13 ,i4, 15,10,17,16,22, 23, 26, 32^ 

In 23 , T.P. Liu studied the following 2x2 simple model which captures the basic 
features of this physical process: 



The numerical issue in such a system is to handle properly the relaxing process 
in its infinite strength limit e — > on coarse computational grids, see 27 - 12 . This 
requirement is twofold as it includes some numerical stability properties and a con- 
sistency with the expected asymptotic behavior in the context of weak solutions. In 
order to derive numerical schemes meeting these criteria, attention has been driven 
onto several kinds of efficient treatments for the relaxation term: see i' 8 * 9 . 15 ' 25 . In- 
deed, robust and efficient approximations can be derived based on "time-splitting" 
algorithms where one alternates at every time-step the resolution of the purely 
convective part and the handling of the reaction process by any O.D.E. solver. 
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x e R,t > with d v f(u,v) < 0. (1.1) 
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It can be sometimes shown rigorously that these numerical schemes converge 
using e.g. the BV compactness framework together with appropriate entropy con- 
sistency properties. Our main goal here is to complete such a study using slightly 
different arguments. Roughly speaking, we are about to consider both the projec- 
tion and the reaction term as being "localized" on a lattice; formally, this amounts 
to going from (1.1) to 



where 5 stands for the Dirac measure in t — 0, At > is a given time-step and 
V h denotes the standard L 2 projection operator on piecewise constant functions 
related to a mesh-size h > 0. Indeed, this new way to handle singular source terms 
in hyperbolic equations can be found in 2 > 7 < 8 > 33 in the context of kinetic equations 
or semilincar relaxation to the one-dimensional conservation law. Such an approach 
based on "measure source terms" allows in particular to recover a fix proposed in 
10 to suppress spurious behavior in the stiff regime. Indeed, the idea proposed in 
10 for reactive chemical flows consists in switching the order in which the reaction 
and projection steps appear inside the Dirac masses, as we shall see in §3.2, §3.3. 

In this paper, we plan to work out only the simpler case of a quasilinear 2x2 
system with relaxation previously considered in 27 > 32 ' 19 j see (2.3). We shall mainly 
rely on an extensive use of some nonconservative products, 5 i 21 i 29 j in order to give 
a rigorous definition of the right-hand side, see Propositions 1 and 2. In return, 
our formulation entails direct stability and convergence proofs exploiting mainly 
estimates easily proved at the level of the original differential problem; this allows 
to bypass some lengthy computations. 

Hence, this work is organized as follows: in §2, we establish some compactness 
properties and study the nonconservative singular system obtained as the relaxation 
and projection terms are concentrated inside Dirac masses, see (2.23). We introduce 
a notion of entropy solution well suited for this kind of problems, see Definition 2. 
Then, in §3, we present some applications to two types of time-splitting numerical 
schemes for which we prove convergence towards the "equilibrium equation" as the 
relaxation parameter diverges, see Theorems 4 and 5. 

2. Study of a singular relaxing system 

2.1. A quasilinear 2x2 model 

We are interested in the Cauchy problem for the following balance laws: 



d t u + d x f{u,v) = (V h {u)-u)S(t-nAt), 




nGN, 



(1.2) 



d t {u + v) + d x f{u) = 0, 

d t v = u, (A(u) - v) , 
u(.,0) = u o ,v(.,0) = v . 



x € R,i > 



(2.3) 
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This system arises e.g. in chromatography, n > 26 > 30 . I n this context, u,v £ [0,1] 
stand for some species densities, one contained in a fluid flowing through a fixed 
bed and the other being absorbed by the material on the bed. This adsorption 
process is modeled by the right-hand side fiR{u, v) — n{A(u) — v), fi 3> 1 where \j [i 
is the relaxation time. This model has been extensively studied within the theory 
of BV functions, 6 ' 34 , by the authors of 31,32 (see also 24 > 19 > 26 > 27 ) under the following 
hypotheses: 

f /(0) = 0, /'(«)>(); 

\ A{0) =0,A(1) = l,A'(u) > 0. 

We recall here some of their main results. 

Definition 1 We say that the pair (u,v) is an entropy solution of (2.3) if: 



(2.4) 



(i) u(.,t) and v(.,t) lie in L 1 n BV(M) for any t £ M + and satisfy (2.3) in the 
sense of distributions; 

(ii) there exists M £ R + such that for any t > 0, s > 0: 

\\u(.,t) - u(.,s)||z,i(R) + \\v(.,t) - v(.,s)\\ L i {R) < M\t-s\; 

(Hi) there holds for any k, I £ R and any nonnegative test-function in 2?(R x R+ ): 



d. 



k\ 



l\ 



/(«) - /(*) 



< fiR(u, v) sgn(v — I) — sgn(u — k) 



(2.5) 



The relevance of this notion comes from the forthcoming theorems whose de- 
tailed proofs are to be found in 32 . We stress however that a more general definition 
of entropy solution to (2.3) has been proposed within the class of L°° functions, 26 . 
Theorem 1 Let u 0} v £ L 1 Pi BV(R); under the hypotheses (2.4), there exists a 
unique entropy solution to (2.3) which satisfies furthermore for any t £ R + : 

(i) /i||i2(u,t;)(.,t)|| I ,i ( R ) < M for some M £ R+; 

(ii) Given another set of initial data uo,wo in L 1 n BV(R), there holds: 

\\u(.,t)-u(.,t)\\ L i(M) + \\v(.,t)-v(.,t)\\ L i(M) < |K-Wo|Ui(R) + |N-6o||l,i(R)- 

The theory of Kruzkov, 18 , ensures that there exists a unique entropy solution 
w to the following "equilibrium" equation since A' > 0: 



d t {w + A{w))+d x f{w)=0, 
w(.,0) = w Q £ L 1 n 



x £ R,t > 



(2.6) 



Theorem 2 Under the assumptions of Theorem 1 and if wq = uq, there exists a 
constant M' £ R+ such that: 



\\u(.,t)-w(. 



M' 

(1? 
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2.2. Introduction of a "localized" singular system 

In all the sequel, At and h will stand for fixed positive parameters related to 
some cartesian computational discretization. 

• We first introduce a projector onto piecewise constant functions; for any j e Z, 
we denote Cj d = [(j — \)h, (j + \)h[. Then, we define: 



V h : L l C\BV(R) -> L l C\ 



(p <-> — / ip(x).dx 



(2.7) 



def 

Let I n = [nAt, (n + l)At[; we define two functions which are Lipschitz con- 
tinuous with respect to time for e > small enough: 



3t»a<(t)={ nAt+ (t^nAt ]Atj te[nAt _ £jnAt[ 



and: 



(n - l)At, te[(n- l)At, nAt - e[ 

(2.8) 



' t - nAt" 

nAt+ [ ) At, t E[nAt,nAt + e[ 



3 t ^ b c (t) = { \ e 

(n + l)At, t e [nAt + e, (n+l)At[ 

(2.9) 

We can therefore consider the following system for any e > 0, u, v in M + within 
the theory proposed in 32 : 

d t u e + d x f(u e ) = u(V h (u e ) - u e )d t a e - uR(u e , v £ )<9 t 6 £ , 

d t v e = v{V h {v € )-v € )d t a e + uR(u € ,v € )d t b e , iGl,i>0 (2.10) 

w e (.,0) = u ,v e (.,0) = va- 
in the limit e — > 0, both a e and & c converge almost everywhere towards the 
piecewise constant function t i-> Af^ n£N F(t — nAt), where Y stands for the 
classical Heaviside distribution. 

We give at once a relaxation estimate for (2.10) which is uniform in e. 

Lemma 1 Under the assumptions of Theorem 1, the entropy solution of (2.10) 
satisfies the relaxation estimate u||i?(u e , u e )c> t £> e 1 1 (iRx K. + ) — ^ f or some C € K + - 

Proof. We follow the classical ideas from 17 > 26 ' 32 introducing a special 
entropy /entropy-flux pair for (2.10): 
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We observe that under the hypotheses (2.4), the function H is strictly convex. 
Therefore, for any positive test-function, the following inequality holds: 

d t r,{u t ,v' i ) +d x q(u € ) < -ii(A(u £ ) -v € )(u £ - A' 1 ^))^ 
+u([V h {u') - u £ )u e + (V h (v e ) - v^A-^v^dta'. 

We integrate (2.11) on M x [0, T] for any T e R+. Using Jensen's inequality, we can 
get rid of the terms involving the projection operator since d t a e > 0: 

f (P h (u e ) - u e )u e d t a e .dx.dt = h [ (V h (u e f - V h ((u e ) 2 ))d t a e .dt < 0. 
JRx[0,T] J[o,t] ^ ' 

We observe now that since A e C 1 (0, 1) is increasing according to (2.4), the deriva- 
tive of its inverse mapping satisfies the following bound: 

H^-Yllco > if-ip- d = f C. 

\\ A \\c° 

And we have for any w, v in [0, 1]: 

(A- 1 (w)-A- 1 (v))(w-v) > C(w-v) 2 . 
Hence the integral involving V h (v e ) is treated the following way: 

f (V h (v e )-v e )A- 1 (v e )d t a e .dx.dt = 

JRx [0,T] 

- [ (V h (v e ) -v c )(A- 1 (V h (v t )) - A- 1 (v e ))d t a e .dx.dt < 

JRx[0,T] 

-C [ (V h (v e ) - v e ) 2 d t a e .dx.dt < 0. 

JRx[0,T] 

Finally, choosing w = A(u), we have for C > 0, the lower bound of (A' 1 )': 

(A(u e ) - v e )(u e - A- 1 ^)) > C(A(u e ) - v e ) 2 . 
Thus we get the expected relaxation estimate (see also 17 ): 

Cii [ (A(u £ )-v £ ) 2 d t b e .dx.dt< [ (r?(u £ , v e )(x, 0) - r?(u £ , v e )) (x, T).dx. 

JRx[0,T] JR ^ ' 

Since BV(R) C L°°(R), we are done. □ 

Another important feature of (2.10) lies in the following compactness result 
which will be of constant use in the sequel of the paper. 

Lemma 2 Under the assumptions of Theorem 1, let (u e ,v e ) be a sequence of en- 
tropy solutions to (2.10). Then (u e ,v e ) is relatively compact in Lj oc (M. x IR+) as 
e -> 0. 

Proof. It follows the classical BV stability canvas and makes use of the 
so-called "quasi-monotonicity" property of the relaxation term, 1 > 26 > 31 . 
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• L 1 (M) bound: we multiply (2.10) by (sgn(u e ), sgn(u e )) T and integrate on x E 
R. The point to be checked is the behavior of the source terms, but we see 
that: 

f sgn{u £ )(V h (u £ )-u e )(x,t).dx<J2 I \u € (x,t)\.dx-\\u%,t)\\ L1{m <0. 
Jr jeZ JCj 

Also, we notice that R(0, 0) = and we can use the mean-value theorem. 
Thanks to the sign assumption on A' in (2.4), we get: 

/ (sgn(i; £ ) - sgn(u e ))(A'(£)u e - v e ).dx < 0. 
Jr 

Since d t a e and d t b e arc nonncgative, we derive finally: 

h e {;t)\\LHR) + \\v e {.,t)\\LHR) < \\u \\ L ^R) + IMUi(R)- (2-12) 



• BV(K) bound: Relying on 18 , we can differentiate the system (2.10) with 
respect to x, we multiply by (sgn(d x u e ), sgn(d x v' i )) T and integrate on i 6 I. 
Using exactly the same arguments, we derive: 

TV(u e )(.,t)+TV{v e ){.,t) <TV(u ) + TV(v ). (2.13) 



• Let t > and t = t + £, £ > 0; we compute 



/ \u e (x,t) - u e (x,t )\.dx < f R f*\d x f(u e )\.dt.dx 
Jr 

+ J R Jt u\R(u*,v*)\d t V.dt.dx 
+ f R J to v\V h (u £ )-u*\d t a*.dt.dx, 

/ \v e (x,t) -v"{x,t )\.dx < f R f* ulRtu^v^ldtb'.dt.dx 
Jr 

+ f R flv\VHv t )-v i \d t a\dt.dx 1 

in conjunction with the bounds: 

iRiu^v^dtb'.dx.dt < 0(1), 

f OV) - « £ | + PV) - v e \)(x,t).dx < h(TV(u ) + TV(v )). 
k Jr v 1 

And we derive: 

f \u'(x,t + £)-u*(x,t)\ \v'(x,t + Z)-v'(x,t)\ 
sup / 1 .dx.dt < 

C#0JRx[0,T] S £ 

(T + At).{(uh + Li P (f))(TV(u a ) + TV(vo)) + 0(1)} 

(2.14) 



Rx[0,T] 
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Therefore, we see that the sequence (u e ,v e ) lies in BV; oc (R x R+) and it remains 
to invoke Helly's compactness principle to conclude the proof. □ 

We notice that because of the Dirac masses in time arising at the right-hand 
side of the equations, layers are likely to appear every time nAt, neN, and 

u £ ,v f - ^C°(R + ;L 1 (R)), e -> 0, 

as in (for instance) 28 in the case of initial data being not "well-prepared" ; or it can 
be evidenced just by taking u = 0, v (x) = sm(wx/h) inside (2.23). 

2.3. A meaning for the ambiguous products 

The preceding bounds for u e imply some compactness for the terms lying in the 
right-hand side in the weak-* topology of measures on R x R+. More precisely, we 
plan to shed some light on the products emanating in (2.10) in the limit e — > 0. 
First, we state a stabilization result. 

Lemma 3 Under the assumptions of Theorem 1, the sequences of entropy solutions 
(u^,v^) to the Cauchy problem with initial data 

(4,«o) e Ll n BV(R) for: 

d t u< + C.0*/(« c ) - uAt{V h {u<) - u<), d t u< + (.d x f(u<) = -nAtR(u<,v<), 
d t vi = vAt(T h (v<) - v<), "' d t vi = ^AtR{u^, v<), 

are relatively compact in Lj oc (R x R+) and belong to C°(M+; L 1 (M)) as ( — > 0. 

Proof. It follows from the same arguments than the proof of Lemma 2. □ 
Proposition 1 Under the hypotheses of Lemma 2, (u e ,v e ) — > (u,v) in L i 1 oc (KxR+) 
as e — > up to a subsequence. Moreover, there holds: 

V h (u C ) — U e \ „ f weak-* M 

1 a f a 



V h (v e ) - v e J * 

|> Of (£$:S )<*'>*)«<'-»"> 

where u,v satisfy the following differential equations for r £ [0, 1]: 

it '.)-•>*( $$Z',). <«<» 



v 



together with the initial data, ieK,i = nAt for neN,: 

u(x,t = 0) = u(x,t - 0), v(x,t = 0) = u(x,t-0). (2.17) 

Proof. We want to compute the value of the following expression for every 
</> £ C"(R x R+), the space of continuous functions with compact support in R x R+: 

f {V h {u e ) - u e )d t a e <j)(x, t).dx.dt = 

r r nAt At ( 2 - 18 ) 

V / / {T h {u e ) -u e )— 4>{x,t).dx.dt. 
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We pick aneN, and define t = 1 + - " At . Thus, according to this "inner variable" , 
the system (2.10) rewrites for r G [0, 1]: 

j d T u e + ed x f(u £ ) = uAt(V h (u e ) -u e ), 
\ d T v £ = vAt{T h {v e ) - v £ ), 

together with the initial data: 

u e (x, t = 0) = u e (x, nAt — e), v e (x, r = 0) = v e (x, nAt — e). 
We perform the change of variables tH-rin (2.18) to get: 

f f nAt At 

(V h (u e )-u e )—<f>{x,t).dx.dt = 



R JnAt-e 



E / / ( :P ' l ( se ) " - " £ )( X ' T ) — <p(x,nAt + 6T-e).dx.edT. 

Lemma 3 implies that u e satisfies (2.15), (2.16) in the limit e — > and we are done. 
□ 

Remark 1 Thanks to the linear form of the projection term, it is possible to define 
directly the limit of (V h (u e ) — u e )dta e as e — > as a distribution of order zero (i.e. 
a bounded measure) provided At = 0(h). Let (p G 2?(K x we can compute: 

J2 f (V h (u)-u)(j){x,nAt).dx= J2 h{v h {u)V h {4>)-V h (u(t>)Y]h,nAt). 

And we have £ j£Z \V h {u)V h (<f>) - V h (u<t>)\(jh,t) < ||^|| c oTV(«)(.,t). 

Concerning the relaxation term, the situation could be different because of its 
nonlinear structure. Nevertheless, we have the following result. 
Proposition 2 Under the hypotheses of Lemma 2, (u e , v e ) — > (u, v) in L ; 1 oc (RxM+) 
as e — > up to a subsequence. Moreover, there holds: 

(A(u e )-v e )d t b eWeak -^* M V At( [ (A(u)-v)(x,T)dAs(t-nAt) (2.19) 
where u,v satisfy the following differential equations for r G [0, 1]: 

<"»> 

together with the initial data, = nAt for neN,: 

u(x,T = 0)=u(x,t-0), v(x,t = 0) = v(x,t -0). (2.21) 

Proof. Let 4> G C°(M x we intend now to pass to the limit e — » in the 
following expression: 

f (A(u e ) - v e )<9 t 6 e <H^ t).dx.dt = 

JRxR+ 

r rnAt+e a . (2.22) 

(A(u e )-v e )—<j)(x,t).dx.dt. 
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We pick a n £ N» and define r = *~" At as an "inner variable" inside the system 
(2.10) which reads for r £ [0, 1]: 

d T u e + ed x f{u e ) = -nAt(A(u e ) - v e ) , 
d T v e =nM(A(u e ) -r), 



together with the initial data: 

u e (x,r = 0) = u € (x,nAt), 



v e (x,r = 0) = v e (x,nAt). 



In order to conclude the proof, we perform the change of variables t i->- t in (2.22) 
and invoke Lemma 3 to derive (2.19), (2.20) in the limit e — > 0. □ 

2.4. A contraction result "a la Kruzkov" 

We can now consider weak solutions of the following singular relaxing system 
for n,i> in R+, x £ R,t > 0: 



d t u + d x f(u) = J2 neNt At u(V h (u) - u) - nR{u, v) 6(t - nAt), 



d tv - E„ eN , At [v{V h {v) -v) + {iR{u, v)J 5(t - nAt), 
u{., 0) = it £ L 1 n BVQBL), v(., 0)=v eL 1 n BV(R), 



(2.23) 



where all the measure- valued source terms are always to be understood in the sense 
of distributions according to the results stated in Propositions 1 and 2. We first 
extend the notion of "entropy solution" to (2.23). 

Definition 2 We say that the pair (u,v) is an entropy solution of (2.23) if: 

(i) u and v He m L 1 n BV(R x [0,T]) for any T £ M+ and satisfy (2.23) m the 
sense of distributions; 

(ii) there holds for any k, I £ M and any nonnegative test-function in 2?(R x M.f): 



\u-k\ + \v- l\ 



/(«) - /(*) 



< 



At v(V h {u) - u)sgn{u - k) + v(P h (v) - v)sgn(v - l)+ ^ 2 ^ 



[iR{u 1 v)(sgn(v — I) — sgn(u — k)) S(t — nAt). 



The main result of this section states that system (2.23), (2.24) is i 1 (M)-contractive. 

Theorem 3 Let (u ,v ) £ L 1 n BV(R); under the hypotheses (2.4), there exists a 
unique entropy solution (u,v) to (2.23). Moreover, if(uo,vo) stands for another set 
of initial data in L 1 D BV(R) and (u, v) for its associated entropy solution, there 
holds for any t £ E+ : 



u(.,t) - u(.,t)\\ L i(m + \\v(.,t) - v{.,t)\\ L i(m < \\u - u \\ L i(m + \\v - v \\ L i {R) . 
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Despite the fact that the solution of (2.23), (2.24) suffer discontinuities on the 
lines t = nAt, n e N», its BV regularity allows one to define a set Af C [0, T] of 
measure zero containing and T such that each point in the strip (x,t) G R x [0, T], 
t £ Af is either an approximate jump or an approximate continuity location, 6 . 

Proof. We need first to check the regularity criterion (i) in Definition 2. This 
is indeed a direct consequence of the uniform bounds (2.12), (2.13), (2.14) in the 
proof of Lemma 2. We start from the entropy inequality (2.5) satisfied by u e , v e for 
e > 0: 



d t 



k\ 



l\ 



/K) - /(*) 



< 



v[(V h (u e ) - u e )sgn(u e -k) + (V h (v e ) - v e )sgn(v e - l)]d t a e 
LiR(u e , v e )(sgn(v e -I)- sgn(u £ - k))d t b e . 



(2.25) 



Therefore, the inequality (2.24) comes in the limit of (2.25) by passing to the limit 
e — > as in the proofs of Propositions 1 and 2. 

Concerning the contraction property which implies uniqueness, we follow the 
ideas of 18 ' 20 > 32 . Let us now introduce a nonnegative test function ty(x,t,y,s) = 
ip(x, t).((x—y).((t—s) belonging to 2?((RxK+) 2 ) with £ the standard approximation 
of the Dirac mass. We test the preceding entropy inequalities (2.25) for both u e ,v e 
and u e ,v e with this particular function, we add and integrate on (R x R+) 2 . For 
e > 0, the standard theory of Kruzkov allows to let ( concentrate to the Dirac 
measure since a e ,b e are Lipschitz functions. We obtain: 



f(u e )-f(u e ) 



< 



v[{{V h {u € ) - u e ) - {V h (u e ) - w £ ))sgn(u £ - u e )+ 
{(V h (v e ) - w £ ) - (V h (v e ) - £ £ ))sgn(t; £ - d e )]d t a e + 
Li{R(u e ,v e ) - R{u\i e )) (sgn{v c - tf ) - sgn(u £ - u e ))d t b e . 



(2.26) 



We study mainly the zero-order terms. Using the mean-value theorem for the 
relaxation term leads to 

/ (sgn(w £ - £ £ ) - sgn(u £ - u e )) (^(u £ ) - A(u e ) - v e + i e \d t b e ^{x, t).dx.dt 
= [ A'(<e)(sgn(« £ -d e )(u e -u £ ) - |u £ - u £ \)d t b e ip{x,t).dx.dt 
+ f (sgn(u £ -u e )(v e -tf) - \v e -v e \)d t b e ip{x,t).dx.dt, 

for a £ e [min(u £ , u £ ), max(w £ , u e )]. Since d t b e > and A' > 0, this term gives 
a nonpositive contribution on the right-hand side of the entropy inequality (2.26). 
Concerning the projection operator, we see that 

/ sgn(u £ - u e )(v h {u e - u £ ) - (u e - u e ))d t a e ^(x 7 t).dx.dt = 

f (sgn(u £ - u e ) - sgn('P' l (u £ - u e ))) {v h (u e - u c ) - (u e - u e ))d t a e ^p.dx.dt 

+ f sgn(V h {u e - u e )) (v h {u e - u e ) - (u £ - u £ )) d t a e ^{x, t).dx.dt. 
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Since V h (u e — u £ ) is constant on each cell Cj, the last integral vanishes. Once again, 
since d t a e > 0, this term gives a nonpositive contribution in (2.26). We end up with: 



u"-u e + v e 



/(« £ ) - f(u e ) 



< 0. 



It remains now to select ip as in 18 ' 20 ' 3 to complete the proof of Theorem 3. □ 
Remark 2 We can notice that this contraction property gives back the BV(M) reg- 
ularity for u(.,t), v(.,t) as soon as u 0} v belong to BV(M.) since the problem (2.23) 
is translation invariant. 

We notice finally that since V h is linear, any weak solution of (2.23) also satisfies: 

d t (u + v) + d x f(u) = EneN, v£t[v h (u + v)-(u + v)] S(t - nAt), 
Q tv = Enefh At ["(P fc («) ~v) + ^R{u 7 «)] S(t - nAt). 

3. Study of the resulting "time splitting" numerical schemes 

3.1. A technical lemma 

We point out quickly that the special structure of the "singular projection term" 
allows us to refine a bit the entropy inequality (2.24). 

Lemma 4 In the sense of Proposition 1, the following holds for any k £ R and any 
nonnegative continuous compactly supported test-function: 



Y vAt(v h (u) - u^j sgn(u - k).5(t - nAt) < 

nGN» 

(\P h {u) - k\ - \u - fc|) (1 - exp(-i/At)) .8{t - nAt). 



Proof. Let < (j> £ C°(R x R+); we perform the following splitting: 

J sgn(u — k) = sgn(P h (u) — k) + (sgn(w — k) — sgn(V h (u) — k)), 
\ V h (u) — u = (P h (u) -k)-(u- k). 

Since we have sgn(x).y < \y\ for any x, y in R 2 , we get out of (2.16) for any k £ R: 

d T \u-k\ <vAt(\V h {u)-k\-\u-k\}. 
We observe that d T \V h (u) — k\ = along the flow of (2.16) and this implies: 

\u-k\{x,l) - \u-k\(x,0) < (1 -cxv(-vAt)).(\V h (u) -k\ - \u - fc|)(x,0). 
Therefore, we just have to recall the following from (2.15) to conclude the proof: 

Y vAi J (^J s S n (" - k )(V h (u) - u){x,r).dT S j <j>(x,nAt).dx = 
^ / ( / d T \u — k\(x, r).dr J <j>(x, nAt).dx. 
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□ 



This means in particular that the following inequality holds as a consequence of 
(2.24) and this last lemma: 



d t 



\u — k\ + \v — . 



/(«)-/(*) 



< 



[(\V h (u) -k\- \u -k\ + \V h (v) -l\-\ v -l\)(l- exp(-^Ai)) (3 2?) 



n£N, 



+[iAtR(u, v) (sgn(v — I) — sgn(w — fc)) <5(t — nAt), 



in the sense of Proposition 2. One can notice the similarity between this last inequal- 
ity and the one derived in 31 §3; indeed the numerical approximation considered 
by these authors matches ours in the special case where v — > +oo and a backward 
Euler solver is chosen instead of (2.20). This reduces to a time-splitting scheme with 
some kind of "straight lines" approximation of the quantity (2.19). In contrast, we 
recall that it has been proved in 3 §5, that one has: 

Y vAt{v h {u) - uj (sgn(?i - k) - sgn(T h (u) - fc)).5(i - nAt) < 0, 



for any positive continuous test-function in the case where u G 



+• r.i 



3.2. The classical time- splitting scheme 

We plan to study the behavior of the system (2.23) in the limit v — > +oo. 
Roughly speaking, this reduces to insert a projection stage every time nAt, neN, 
before igniting the relaxation mechanism. 

Proposition 3 Under the assumptions of Theorem 3, the sequence of entropy solu- 
tions of (2.23) converges strongly in Lj oc (R x M+) as v — > +oo towards the entropy 
solution of: 



d t u + d x f(u) = - J2neN, pAtR(u, v).S(t - nAt), 
®tv = J2 n eN, V&tR(u, v).S(t - nAt), 



(3.28) 



where the right-hand side is still defined by means of (2.19), (2.20), but with the 
piecewise constant initial data for t = nAt, neN,: 



u(x, t = 0) = V h {u){x, t-0), v(x, t = 0) = V h (v)(x, t-0). 



(3.29) 



Moreover, the following entropy inequality holds true for any nonnegative test func- 
tion in V(R x 



d t [\u - k\ + \v - l\ 

EneN t {(\-P h ^)-k\-\u- 



f(u)-.f(k)\ < 



k\) + (\V h (v)-l\ - \v-l\) + 



(3.30) 



fiAt.R(u,v)(sgn(v — I) — sgn(u — fc))|<5(t — nAt), 
with k,l inR and according to (2.19), (2.20), (2.21). 
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Proof. We split it into several steps for the sake of clarity, 
(i) Construction of the solution of (2.16) by a fixed point argument; for r G [0, 1], 
we consider any function w G C°(0, 1; together with: 

d T u = i>At{T h {w) - u), w(.,0)GL 1 (M) 

whose solution is obvious: 

u(x, t) = exp(— vAtr)u(x, 0) + vAt cxp(— vAtr) I ex${vAts)V h {w)(x, s).ds. 

Jo 

This defines a linear mapping w H» u in C°(0, 1; i 1 (R)) for any choice of the initial 
datum. We want to establish a contraction property: for all Wi,w 2 and any r G 
[0, 1], we observe that 

f \ui ~ u 2 \(T).dx < vAt f exp(vAt(s - r)) f \V h (w 1 ) - P h (w 2 )\(x, s).dx.ds 
Js, Jo Jr 

< (1 - exp(-i/Ai))||u;i - w 2 \\c»(os;L^m)) 

By Picard's fixed point theorem, we obtain a unique solution in C°(0, 1; i 1 (R)) for 
any v G R + to: 



d T u = vAt{V h {u)~u) 1 u(.,0)eL 1 n 

Now, since d T V h {u) — along its flow, this differential equation therefore admits 
an explicit solution: 

u(., t) = exp(— vAtr)u{., 0) + (1 - exp(-uAtT))V h (u)(., 0). 
In particular, we get: 

(u - V h {u^j (., r) = exp(-^Atr) (u - V h {u^j (., 0). (3.31) 
(ii) Limit as v — >• +oo of the measure (2.15); for any <j> G C°(R x R+), we have: 

"At [ ([ (V h {u) - u){x,t)At\ cf>(x,nAt).dx = 
vAt j (T h (u)-u)(x,nAt-0) ( / exp(-i*-).dr ) (j>{ Xl nAt).dx 



u ^4°° J2 I (V h (u) - u)(x,nAt - 0)cj)(x,nAt).dx, 



which is a finite sum provided At remains strictly positive. 

(iii) Strong compactness as v — > +oo: the bound (2.14) seems to blow up despite 
the fact that (2.12), (2.13) still hold. Making use of (3.31) in the last step of the 
proof of Lemma 2, we notice that for any T G R + and £ > 0: 

1+[T/At] 



/ V vAtl / u-V h {u) + v-P h (v) (x,T)d T \ .dx.ds < 

Js.x[t,t+£] n=1 [Jo J 

^(T + At)(TV(«o) + TV(v j) (1 - exp(-i/At)) . 



(3.32) 
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Therefore, exploiting (2.15) and Lemma 1, we can deduce a fine estimate for the 
weak solutions of (2.23), (2.15), (2.16) which is uniform in both v and [i: 

sup / l"0M + Q -*(*,*)! + KM + fl-«(*,t)l -(fadt < 

?#oJkx[o,t] ? ? (3.33) 

(t + At){(VAt + l»p(/))(tv(uo) + ry(« )) + 0(1)} 

Together with (2.12), (2.13), this provides compactness in Lj oc (R x M+) by means of 
Helly's theorem. And in the limit, we get the following system: (recall from (1.2)) 

d t u + d x f{u) = Y^nen, - «)(■, nAt - 0) - »AtR(u, vj\ 5(t - nAt), 

d tv = E„ GN , [CP h ( v ) - ")(-. ™ A< - °) + MAti?(u, w)] <J(t - nAt). 

Since we know that a £ (t) — b € (t — e), we must handle the projection step first and 
the initial data for (2.20) becomes (3.29). The preceding system then rewrites like 
(3.28). 

(iv) Entropy inequality: Using the same arguments together with Lemma 4, we 
derive for any fcel: 

V vAt [ ([ sgn(u - k)(V h (u) - u)(x,T).d,T) <j>(x,nAt).dx = 

^ / (\u-k\(x,l)-\u-k\(x,0)^<f>(x,nAt).dx < 

nes, ^ R 

X! / (\'P h (u)-k\-\u-k\Yx,nAt-0)(l-cx-p(-i'At))(l)(x,nAt).dx 
""4°° I (l^iu) - k\-\u- k\y x ,nAt-0)(t>(x,nAt).dx. 

n£K, ^ R 

This way, we get (3.30) out of (2.24) and the contraction property of Theorem 3 is 
also preserved. □ 

One can observe that as soon as the initial data of (3.28), (2.19), (2.20), (3.29) 
is ^-invariant and the CFL condition Ai||/'(u)||z,oo < h holds, its entropy solu- 
tion matches the piecewise constant approximation generated by a classical time- 
splitting Godunov scheme, see e.g. 1,31 and 9 . 

Theorem 4 Under the assumptions (2.4), the CFL restriction Ai||/'(ti)||.i,°° = h 
and for v = A(u ), u e L 1 n BV(M), the entropy solutions of (3.28), (3.30), 
(2.19), (2.20), (3.29) converge strongly in Lj oc (R x R+) towards u, the entropy 
solution in the sense of Kruzkov to: 

d t (u + A(u))+d x f(u) = 0, u(x,0) = u . 

as fiAt — > +oo, h — > 0. 

Proof. As we did before concerning Proposition 3, we are about to split the 
forthcoming proof for the sake of clarity. 
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< +00 



(i) Study of the solution of (2.20); for t £ [0, 1], we rewrite its second equation 

as: 

d T v = /j,At(A(v) - v), A(v) = A(u(x,0) +v(x,0) -v(x,t)). (3.34) 

We follow the classical approach, see e.g. 4 , and we introduce the following Banach 
space for some k £ M+ to be fixed later: 

X = < ip £ C°(0, 1; i 1 (M)) such that sup ( cxp(-k.t)\\ip(. 7 t)\\ L i( R) ) 
{ te[o,i] v ' 

For any id e I, we consider the mapping lo^iiel, 

v(x, t) = v(x, 0) exp(— ^Air) + / fiAt exp(/iAt(s — t))A(w)(x, s).ds, 

Jo 

which turns out to be a contraction for k > fiAt(Lip(A) — 1). By Picard's fixed 
point theorem, we get a unique solution to (3.34) v £ X for any [i £ M+. Moreover, 
it satisfies 

d T (A(v) - v)(x, t) = -fj,(l + A'(u{x, 0) + v{x, 0) - v(x, t))) (A(v) - v)(x, r), 
which leads to: 

\A(u) — v\(x, t) < exp(— /zAtr)| A(u) — v\(x, 0). 

(ii) Behavior of (3.30) as ^ — > +oo; the bound (3.33) can be simplified by means 
of the CFL condition: 

sup / l"(*,* + Q -«(*,*)! + NM + O-^t)! ^ < 

«#oJrx[o,t] £, £ (3.35) 

(T + At){2Li P (f)(TV(uo) + TV(voj) + 0(1)}. 

This provides strong Lj oc compactness and for any positive continuous test function 
4> £ C°(R x R+), we have: 

At fj,( / (A(u) — u)(sgn(w — /) — sgn(u — k))(x,T).dT\(f>{x,nAt).dx 
<2V [ \A(V h (u))-V h (v)\{x,nAt)( [ fiAt. exp(- ^Atr) .dr) (f>(x, nAt).dx 



nEN, 



<o r' cx ir At) )u\\c°. 



fiAt 

We used the preceding estimates in conjunction with the bound in Lemma 1. There- 
fore, this term vanishes in the limit /j,At — > +oo, /i — !• 0. Moreover, by Jensen's in- 
equality for convex functions, we also have thanks to (2.13) and the CFL condition, 

I (\P h {u)-k\-\u-k\^{x,nAt-Q)(l)(x,nAt).dx< 
I \V h (u)-u\(x,nAt-0)\V h (4>)-4>\(x,nAt).dx< 



hEnen, &tLip(f)(TV(u ) + TV(vo))Lip(4>) h 3° 0. 
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(iii) It remains to pass to the limit in (3.30); we perform the following splitting 
for x G M and t G [nAt, (n + l)At[ 

\v - A(u)\(x,t) < \v(x,t) - v(x,nAt)\ + \v - V h (v)\(x,nAt) + 
\A(V h (u))-V h (v)\(x,nAt)+ 

\A(V h (u)) - A{u)\{x,nAt) + \A{u){x,t) - A{u){x,nAt)\, 

each term converging to zero in L 1 as [iAt — > +00, h — > under the prescribed 
CFL condition. Therefore, we fix I = A(k) G R and we observe that for any positive 
4> G X>(R x R+) 

/ |w - l\d t ip(x,t).dx.dt -)■ / — A(fe)|atV(a:,t)-<ia:-d*, 

which leads to the integral form of Kruzkov's entropy inequality: 

/ fc| + - A(/c)|)d t V>0M) + f(u)-f{k)d x tp{x,t).dx.dt>0. 

□ 

3.3. A modified time- splitting scheme 

It is clear that the estimates of Lemmas 2 and 1 still hold for the slightly modified 
system (compare with (2.10)): 

dtu e + d x f(u e ) = -fiR{u*,v e )d t a e + v{V h {vf) - u e )d t b e , 
d t v e = nR(u e ,v e )d t a e + v(V h {v e )-v e )dtb e , xeR,t>0 (3.36) 

u e (.,0) = u o ,v e (.,0) = v Q . 

The Cauchy problem for (3.36) still can be studied within the framework of 32 for 
any e > 0. According to the same ideas, we send e — > to get the following result. 
Proposition 4 Let (uo,Vo) G L 1 n BV(R); under the assumptions (2.4), the se- 
quence (u e 7 v e ) of entropy solutions to (3.36) converges in Lj oc (R x R+) towards the 
entropy solution in the sense of Definition 2 to 



d t u + d x f(u) = J2 n&h At - nR(u, v) + v(V h { 



[u - u 



S(t - nAt), 



(3.37) 



dtv = Enefh At [/ifl(u, v) + v(V h (v) - v) j S(t - nAt) , 
which satisfies for any k, I G R and any nonnegative test function in X>(R x R+): 



d t [\u-k\ + \v-l\\ +d x \f(u)- f(k) 
At /j,R(u,v)(sgn(v — I) — sgn(u — k)) 



(3.38) 



+v{T> h (u) - u)sgn(u - k) + v(V h {v) - v)sgn(v - I) S(t - nAt) 



in the sense of Propositions 1 and 2. Given any other set of initial data {uq,Vq) in 
L 1 n BV(M), this last inequality implies the contraction property for any t G R + : 

\\u(.,t) - u(.,t)|| L i (R) + \\v(.,t) - v(.,t)\\ L i(jn) < \\u - uo||li(r) + \\v - voIIli(r)- 
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Proof. Similar estimates than the ones shown in the proof of Lemma 2 ensure 
the strong L] oc compactness. The inequality (3.38) comes from (2.25) together with 
Propositions 1 and 2. The contraction property in i 1 (R) is proved the same way 
than for Theorem 3. □ 

Remark 3 It is important to notice that the entropy solutions of (3.37) do not 
coincide with the ones of (2.23) since the relaxation and the projection steps are not 
applied in the same order. This means in particular that one has to be careful when 
dealing with the right-hand sides as different processes are concentrated inside the 
Dirac masses. Hence there is no uniqueness problem for both systems (2.23), (3.37) 
relying on Propositions 1 and 2 together with the inequalities (2.24) an d (3.38). 

We plan now to carry out the same program as we did for the classical time- 
splitting scheme in the preceding section. By letting v — > +oo, we derive an equation 
whose entropy solution coincide with the piecewise constant approximation gener- 
ated by a numerical scheme designed in the context of reactive Euler equations for 
which the reaction process operates before the projection stage, 10 . 
Proposition 5 Assume (2.4) and («o, «o) ^ L 1 C\ BV(M), the sequence of entropy 
solutions to (3.37) converges strongly in Lj oc (R x R+) as v — > +oo towards the 
entropy solution of: 

d t u + d x f(u) = J2 neN , [ - nAtR(u, v) + (V h (u) - u)(., nAt + 0)] S(t - nAt), 
d tv = E„ eN , [v&tR(u, v) + (T h (v) - v)(., nAt + 0)] 6(t - nAt), 

(3.39) 

where the right-hand sides are defined by means of (2.19), (2.20), (2.21). More- 
over, the following entropy inequality holds true for any nonnegative test function 
in D(R x R+) with k, I in M: 

d t [\u-k\ + \v-l\\ +d x \f(u)- f(k)\ < 
E„ eN , {»At.R(u,v)(sgn{v -I)- sgn(u - k))+ (3.40) 
(\P h (u) - k\ - \u - k\) + (\V h (v) - l\ - \v - l\)}8(t - nAt), 

The final statement is concerned with the relaxation limit in (3.39). 
Theorem 5 Under the assumptions (2.4), the CFL restriction At\\f'(u) = h 
and for v = A(u ), u £ L 1 D BV(R), the entropy solutions of (3.39), (3.40), 
(2.19), (2.20), (2.21) converge strongly in Lj oc (R x R+) towards u, the entropy 
solution in the sense of Kruzkov to: 

d t (u + A(u))+d x f(u)=0, u(x,0)=u o . 

as [iAt — > +oo, h — > 0. 

The proofs of these last two statements are completely similar to the ones of 
Proposition 3 and Theorem 4 hence we skip them. 



4. Conclusion 
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In this paper, we considered some "measure source terms" for a quasilinear relax- 
ation system whose weak solutions coincide with piecewise constant approximations 
generated by commonly-used time-splitting numerical schemes. This allows to es- 
tablish convergence results directly relying on properties of the underlying original 
system and bypassing some heavy computations required to show the stability of 
the numerical processes (see e.g. §2.5 in 7 ). 
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